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We study cold atomic gases with a contact interaction and confined into one-dimension. Crossing 
the confinement induced resonance the correlation between the bosons increases, and introduces an 
effective range for the interaction potential. Using the mapping onto the sine-Gordon model and 
a Hubbard model in the strongly interacting regime allows us to derive the phase diagram in the 
presence of an optical lattice. We demonstrate the appearance of a phase transition from a Luttinger 
liquid with algebraic correlations into a crystalline phase with a particle on every second lattice site. 



Cold atomic gases confined into one-dimension exhibit 
remarkable properties as the interplay between interac- 
tions and reduced dimensions strongly enhances quan- 
tum fluctuations. The most prominent example is the 
appearance of a Tonks-Girardeau gas for bosonic parti- 
cles [l|, |2j], and the possiblity to pin the bosons into a 
Mott insulating phase for arbitrary weak optical lattices 
3 > 13 ■ Most remarkably, it has recently been proposed 
5] and experimentally observed |6| , that it is possible to 
access a regime, where the bosonic many body system 
exhibits even stronger correlations. This opens the ques- 
tion, whether it is possible to enhance the correlations 
to a point, where the bosonic systems forms a crystalline 
ground state. In this letter, we demonstrate that indeed 
in the presence of an optical lattice a solid phase appears. 

The transverse confinement for cold atomic gases is 
experimentally efficiently achieved using optical lattices 
[H, 0] or atomic chips Q . Within this one-dimensional 
regime with the kinetic energy of the particles much lower 
than the transverse trapping frequency, the interaction 
between the particles is described by the one-dimensional 
scattering length a 1D @ . Remarkably, the system can un- 
dergo a confinement induced resonance, where the scat- 
tering length crosses zero. For a 1D < 0, the properties 
of the system have been studied in terms of the exactly 
solvable Lieb-Liniger model 11], while at a lu = 
the system is denoted as Tonks-Girardeau gas. Cross- 
ing the confinement induced resonance with a 1D > the 
mathematical model describing the system admits a two- 
particle bound state. Then, the physical state smoothly 
connected to the Tonks-Girardeau gas corresponds to an 
highly excited state of the mathematical model; a regime 
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denoted as Super- Tonks-Girardeau gas 

In this letter, we analyze the phase diagram within this 
regime and demonstrate the appearance of a solid phase 
in the presence of an optical lattice with a bosonic par- 
ticle on every second lattice site. A simplified picture of 
this transition is that the particles acquire a finite range 
interaction ~ 

Gid and essentially behave as hard spheres 
Then, it is natural to expect the appearance of a 
solid phase for a density comparable to the range of the 
interaction. The rigorous derivation of the phase diagram 
follows in two steps: First, we analyze whether an arbi- 
trary weak optical lattice allows to pin the solid structure. 
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FIG. 1: (a) Phase diagram: The solid phase appears at in- 
termediate strength of the optical lattice for ai D n > 0.2. The 
blue (dashed) line derives from the transition within the Hub- 
bard model, while the green (dotted) line describes the lower 
bound for the solid phase predicted from the sine-Gordon 
model, (b) Illustration for the two degenerate ground states 
with an atom on every second lattice site, (c) Luttinger pa- 
rameter derived from the exact Bethe ansatz equation. The 
dashed line denotes the asymptotic behavior K = (1 



Using the mapping to the exactly solvable sine-Gordon 
model, we find, that a finite strength of the optical lattice 
is required. Therefore, we focus on deep optical lattices 
in a second step, and provide the derivation of a Hub- 
bard model close to the confinement induced resonance. 
The combination of the two methods allows us to identify 
an experimentally accessible region, where a solid phase 
appears, see Fig. Q] 

We start with the many-body theory describing 
bosonic particles confined into one-dimension. Intro- 
ducing the bosonic field operators tf)'(x) and ip(x), the 
Hamiltonian takes the form 



H n = 



dx ijr (x) 



ti 2 
-—A 
2m 



V(x) 
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dxdy U B (x - y)if> 1 (x)^ (y)ip(y)ip(x). 



Here, V(x) = Vocos 2 (xfc) accounts for the optical lat- 
tice along the tubes. The interaction potential be- 
tween the bosons confined into the lowest state of 
the transverse trapping potential reduces to U B (x) — 
g B S(x) with the coupling strength g B — — 2?i 2 /( ma 1D ) 
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[9(. Here, the one-dimensional scattering length a 1D = 
—a\/a s (1 — Ca s /a±) is related to the three-dimensional 
s-wave scatterling length a s and the transverse confin- 
ing length a_L with C « 1.46 [§]. The system exhibits 
a confinement induced resonance at a s = a±/C, where 
the coupling strength diverges and eventually changes its 
character from repulsive to attractive. 

A physical interpretation of the confinement induced 
resonances is provided by the following property: The 
ID scattering length a 1D describes the distance, where 
the scattering wave function for two particle crosses zero. 
While for a lu < 0, the zero appears in the unphysical re- 
gion < 0, the scattering wave function exhibits a node 
for a 1D > 0. This behavior is achieved by an attractive 
interaction potential U B (x) giving rise to a bound state. 
Then, the scattering wave function is orthogonal to the 
bound state and consequently exhibits a node. Note, 
that the sudden appearance of a bound state is an ar- 
tifact of the mathematical model Eq. ([T]), which is only 
valid in low energy sector with the relevant momenta q 
satisfying the condition qa± < 1. In the physical system 
a bound state is always present; its position across the 
confinement induced resonances has been studied in de- 
tail [Hj]. As a consequence, the atomic system is for all 
values of a 1D a highly excited metastable state, and losses 
via three-body recombination reduce the life time of the 
atomic gas. This indicates that the transition from the 
regime with repulsive interaction into the Super- Tonks- 
Girardeau gas is described by a smooth cross-over. 

In the following, we first focus on the limit of a very 
weak optical lattice Vq <C E r . Then, the low energy prop- 
erties of the strongly interacting bosonic system are well 
described within the hydrodynamics description with the 
bosonic field operator %p(x) ~ y/n + d x 9/n expressed in 
terms of the long-wavelength density and phase fields 
9(x) and 4>(x). The fields satisfy the standard commu- 
tation relation [d x 9(x), (f>(y)] = in5(x — y). The effective 
Hamiltonian in absence of an optical lattice reduces to 



Hn = 



hv. 



dx 



(2) 



The dimensionless Luttinger parameter in the strongly 
interacting regime j B = g B m/nh 2 >• 1 reduces to K = 
(1 — na 1D ) 2 |l0|. This expression remains valid in the 
strongly repulsive situation with a s < 0, as well as in 
the attractive case a s > for |na 1D | <§C 1 [1J, |l5[. In 
the latter case, the dimensionless parameter K < 1 re- 
duces below the non-interacting Fermi limit (K = 1). 
Usually this regime can only be reached for bosonic par- 
ticles through an interaction potential with a finite range. 
Here, such a finite range is achieved from the potential 
U B (x) by the presence of a bound state and the asso- 
ciated node in the two-particle scattering wave function. 
The behavior of the Luttinger parameter K for larger ID 
scattering lengths can be derived from the exact Bethe 
Ansatz equation 15j, see Fig. [TJ Note, that this behav- 



ior is in disagreement with variational Monte-Carlo sim- 
ulations predicting an instability for na 1D > 0.38 and a 
minimal value K k, 0.85 [5j. 

Within this hydrodynamic description the weak opti- 
cal lattice is a relevant perturbation at commensurate 
fillings. Here, we are interested in densities n — I /(ma) 
with a = n/k the lattice spacing and m € N an integer. 
Then the Hamiltonian accounting for the optical lattice 
Vq cos(fcx) takes the form 0, EH 



dx cos (2m6) 



(3) 



with u — KVq/ E r (a/2a) 2 and a a short distance cut- 
off (the cut-off is in the range of the interparticle dis- 
tance a ~ !/"■)■ The low energy description of the in- 
teracting bosonic system H ctt = Hq + H lMicc reduces to 
the quantum sine-Gordon model. This model is exactly 
solvable and exhibits a quantum phase transition from 
a gapless phase with algebraic decay in the superfhud 
correlation function (ifl (x)ip(0)) ~ x~ 1 / 2K as well as in 
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the solid correlation (n(x)n(0)) ~ cos(27rna;)/a; 
gapped and incompressible insulator with long range or- 
der (n(x)n(O)) — n 2 ~ cos(27rn:r). Below the critical 
value K < K m = 2/m 2 , the transition appears for ar- 
bitrary strength of the lattice potential, while for a fixed 
value of u, the transition appears at the universal value 
K = 2/m 2 . Here, K denotes the renormalized Luttinger 
parameter due to the optical lattice; for weak optical 
lattices it is related to the microscopic value K via the 
Kosterlitz-Thouless renormalization group flow (see [TtJ 
for a review). For a bosonic density equal to the lattice 
spacing, i.e., n = k/ir with m = 1, the phase transition 
takes place from the superfluid to the Mott insulating 
phase and has been previously discussed 

In the regime with a positive ID scattering length 
a 1D > 0, it is now possible to access values K < 1. This 
opens the question, whether it is possible to reach the sec- 
ond instability with m — 2 and particle density n = k/2ir, 
i.e., on average there is one bosonic particles distributed 
over two lattice sites. Then, the phase transition takes 
place from a Luttinger liquid with algebraic correlations 
to a crystalline phase. In addition to an excitation gap 
and the incompressibility, the crystalline phase is char- 
acterized by a long range order with a bosonic particle 
localized in every second lattice site. The ground state 
breaks the discrete translation invariance of the system 
and is two-fold degenerate. This property distinguishes 
the solid phase from the Mott insulator at integer fillings. 

The criticial value of the Luttinger parameter, where 
an arbitrary weak optical lattice allows to pin the bosonic 
crystalline structure reduces to K 2 = 1/2. As dis- 
cussed above, this regime can not be accessed. However, 
the optical lattice increases the correlations between the 
bosonic particles. Using the Kosterlitz-Thouless renor- 
malization group flow to lowest order in u for the tran- 
sition line, i.e., K = (1 + it)/2, we can expect the phase 
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transition into the solid phase for a finite strength of the 
optical lattice, see Fig. [TJ For values of the optical lat- 
tice V ~ E r , the effective low energy theory Eq. @ is 
no longer valid, and different approach is required for 
analyzing the appearance of the solid phase. 

In the regime of a strong optical lattice Vq > E r , the 
suitable approach is to map the bosonic system to a Hub- 
bard model and analyze the transition within this regime. 
In the strongly correlated regime with 7 B 3> 1 the con- 
ventional derivation of the Hubbard model fails. Here, 
we first apply an exact mapping of the strongly interact- 
ing bosonic system onto a weakly interacting Fermi gas; 
in the latter situation, the conventional derivation of the 
Hubbard model is valid and allows us to derive the phase 
diagram for strong optical lattices. 

This duality transformation of the strongly interacting 
bosons onto weakly interacting fermions has been pio- 
neered in the past IH, Ijjj]. On the two particle level, it 
requires that the scattering wave function ip B (x) between 
two bosons with the interaction potential U B , is described 
by the a fermionic scattering wave function ip F (x) with a 
novel interaction potential U F via V'b(^) = sgn(x)^ F (cc) 
(here, x denotes the relative coordinate). This property 
is uniquely determined by the pseudo-potential 

W*\<t>)= lirn f [^(e) + ^'H]*[^(e)+^H] (4) 

with g F = 2h 2 a 1D /m the coupling strength and ip' = d x ip 
(</>' = d x 4>) the derivatives of the wave function. It is im- 
portant to note that the role of the ID scattering length 
a 1D is reversed in fermionic pseudo-potential U F as com- 
pared to the bosonic one U B . As a consequence, this 
mapping allows us to transform a strongly interacting 
bosonic model onto a weakly interacting Fermi system. 
Note, that the lim c _ i .Q+ is required in order to avoid a ul- 
traviolet divergence when applying the interaction poten- 
tial on the Greens function. This behavior is in analogy 
to the well known regulariztion of the pseudo-potential 
for 3D s-wave scattering. 

Extending this two-particle analysis to the many- 
body system, therefore maps the bosonic Hamiltonian 
in Eq.([l|) onto a fermionic model 
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V(x) 



ip F (x) 
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dxdy U F (x - y)ipt(x)ipt(y)ip F {y)^ F (x) 



with the fermionic field operators ipi and ip F (x) . The pa- 
rameter 7 F characterizing the strength of the interaction 
in the fermionic model is given by the ratio between the 
kinetic energy E Un — h 2 n 2 /m and the interaction energy 
E iat = n 3 g F , i.e. 7 F = E irit /E hirl = 2na 1D = -1/7 B - The 
ground state wave function \g F ) of the fermionic problem 
is related to the ground state of the bosonic problem \g B ), 

(xt, . . .,x N \g B ) = A(xx,. . .,x N )(xi,. . .,x N \g F ) (6) 
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FIG. 2: Tunneling amplitude 4J (red) and the Wannier func- 
tion overlap \ (blue) for different strengths of the optical 
lattice. The inset shows the renormalization of the nearest- 
neighbor interaction V Q a for large ID scattering lengths ac- 
counting for the influence of higher bands and the proper 
treatment of pseudo-potential U F the at Vo = 4E r . 



with the total asymmetric factor A(x±, . . . ,xn)- For 
bosons with a lu = 0, this mapping reduces to the well 
known relationship between impenetrable bosons and 
fermions in ID [11]. However, here the situation is gener- 
alized to arbitrary strength of the interaction potential. 

In the interesting regime with strong interactions be- 
tween the bosons |7 B | = |1/7f| S> 1, the fermionic system 
is weakly interacting and the conventional approach to 
derive the Hubbard model is valid. Therefore, we obtain 
for a deep optical lattice the fermionic Hubbard model 



Hum — - J ^ 



eje 



C i C ] C j C i 



(7) 



with the fermionic creation (anihilation) operator ct (cj). 
In addition, the hopping amplitude J accounts for the 
single particle band structure tu — — 2 J cos ka, while the 
fermionic pseudo-potential U F gives rise to a dominant 
nearest-neighbor interaction 



2 a 
V = -^E r - 



Vo 

a V E r 



(8) 



Here, x is determined by the overlap between the Wan- 
nier functions w(x) on neighboring lattice sites, 

X = a 3 Jdx \d x w(x)w(x—a) — w(x)d x w(x— a)\ 2 . 

The hopping amplitude J as well as the dimensionless 
overlap x can be efficiently determined numerically for 
different strengths of the optical lattice, see Fig.[5J Note, 
that additional interaction terms are strongly suppressed 
due to the fast decay of the wannier functions. 

At half filling with one particle on every second lat- 
tice site, the Hubbard model Eq. exhibits a quantum 
phase transition from a phase with algebraic correlations 
between the fermions for J 3> V to a charge density wave 
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with an excitation gap for 7> J. The latter phase cor- 
responds to the interesting crystalline phase. The critical 
point for the phase transition is determined by the spe- 
cial point at J = V/2, where the system becomes SU(2) 
invariant and maps to the spin-1/2 Heisenberg model. 
It is this enhanced symmetry which fixes the transition 
point to J = V/2 even in the one-dimensional situation. 

From the behavior of V and J for different strengths of 
the optical lattice, we can now derive the complete phase 
diagram, see Fig. [TJ for very deep optical lattices the 
nearest neighbor interaction is strongly suppressed com- 
pared to the hopping term, see Fig. [2l and consequently, 
the ground state is determined by a Luttinger liquid 
phase with algebraic correlations. Reducing the strength 
of the optical lattice, the nearest-neighbor interaction in- 
creases and a phase transition into the solid phase takes 
place for sufficiently strong interaction a lri n > 0.2. For 
even weaker optical lattices, the mapping to the Hubbard 
model breaks down, and the effective theory is given by 
the sine-Gordon model. The sine-Gordon model requires 
a finite strength of the optical lattice for the appearance 
of the solid phase. Therefore, a second phase transition 
takes place for decreasing optical lattice, and the system 
enters again the Luttinger liquid phase, i.e., the system 
exhibits a remarkable reentrant feature. Consequently, 
we predict the existence of a solid phase for cold atomic 
gases at strong interactions a 1D n > 0.2 and intermediate 
the optical lattices V ~ 3E r . 

Finally, we have to verify the validity of the Hub- 
bard model in the interesting regime with na 1D > 0.2. 
The derivation of the Hubbard model involves two ap- 
proximations: (i) first, we restrict the analysis onto 
the lowest Bloch band, i.e., we introduce a high energy 
cut-off A > a determined by the lattice spacing, (ii) 
Second, the interaction potential U F is treated without 
the proper regularization. The influence of these two- 
approximation has recently been studied in detail for the 
derivation of the Hubbard model in a three-dimenionsal 
optical lattice [20|. Here, the situation is equivalent and 
the main results can be directly carried over. It fol- 
lows, that the Hubbard model is correct for weak in- 
teractions din -C a, while in the interesting parameter 
range a 1D n ~ 0.2 corrections from higher bands and 
the proper treatment of the interaction potential ap- 
pear. The main influence is a renormalization of the 
nearest neighbor interaction strength, which takes the 
from V eB — V/(l + r]V/E r ). Here, rj describes a dimen- 
sionlcss parameter which in general derives from a full 
numerical analysis. However, due to the duality map- 
ping between the Bosons and Fcrmions, we know that 
in the limit a lu /a — > oo the system has to reproduce the 
scattering of non-interacting bosons. This condition fixes 
the parameter to rj = —E r /2J. Therefore, we find that 
the influence of higher bands and the proper treatment 
of the interaction potential increases the strength of the 
nearest-neighbor interaction, see Fig. [2] Therefore, we 



expect that the solid phase appears even for weaker in- 
teractions than shown in Fig. [T] 

The experimental setup required for the observation 
of the solid phase can be achieved by the combination 
of strong transverse confining by an optical lattice with 
a Feshbach resonance to tune the strength of the s-wave 
scattering length. Such a setup has recently been real- 
ized for the observation of correlations beyond the Tonks- 
Girardeau regime Q . An additional weak optical lattice 
along the tubes then opens the path to the experimental 
search of the solid phase. Finally, it is important to note, 
that the behavior of losses by crossing the confinement 
induces resonance are not yet well understood. However, 
for increasing ID scattering length, additional terms to 
the Hamiltonian breaking the integrability of the model, 
e.g., corrections from higher transverse states and addi- 
tional non-universal three-body interactions, provide a 
decay rate and eventually an instability of the Super- 
Tonks-Girardeau gas towards the formation of bound 
states; such a behavior was observed within the varia- 
tional Monte Carlo simulations Q. This implies a finite 
lifetime for the realization of the experiments and sug- 
gests that the search for the solid phase should be per- 
formed for intermediate interaction strengths na lu ~ 0.4. 
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